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Abstract. In a previous paper, using an equation of state of dense matter representing a fluid of Skyrmions 
we constructed the corresponding non-rotating compact-star models in hydrostatic equilibrium; these are mostly 
fluid stars (the Skyrmion fluid) thus naming them Skyrmion Stars. Here we generalize our previous calculations 
by constructing equilibrium sequences of rotating Skyrmion stars in general relativity using the computer code 
RNS developed by Stergioulas. We calculated their masses and radii to be 0.4 < M/Mq < 3.45, and 13.0 km < 
-R < 23.0 km, respectively (R being the circumferential radius of the star). The period of the maximally rotating 
Skyrmion stars is calculated to be 0.8 ms < P < 2.0 ms. We find that a gap (the height between the star 
surface and the inner stable circular orbit) starts to appear for M ~ 2.OM0. Specifically, the Skyrmion star mass 
range with an existing gap is calculated to be 1.8 < M/Mq < 3.0 with the corresponding orbital frequency 
0.8 kHz < visco < 1.3 kHz. We apply our model to the 4U 1820-30 low mass X-ray binary and suggest a plausible 
Skyrmion star candidate in the 4U 1636-53 system. We discuss the difficulties encountered by our model in the 
4U 0614-1-09 case with the highest known Quasi-Periodic Oscillation frequency of 1329 Hz. A comparative study 
of Skyrmion stars and models of neutron stars based on recent/modern equations of state is also presented. 
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1. Introduction 

A key input in determining the structure of compact ob- 
jects is the equation of state (EOS) of high density mat- 
ter whose determination remains a formidable theoreti- 
cal problem. There is no general agreement still on the 
exact composition of dense matter, and on its EOS, es- 
pecially for densities in excess several times nuclear mat- 
ter density (in a regime where the EOS is poorly con- 
strained by nuclear data and experiments). In Ouyed 
& Butler (1999, hereafter OB), we presented and dis- 
cussed the Skyrme model (Skyrme 1962a&b) and its plau- 
sible link to Quantum Chromo-Dynamics (QCD). We ex- 
plained how the connections are intriguing enough that 
the model has been revived in order to study its pre- 
dictions of hadronic interactions. Following an approach 
outlined in Kalbermann (1997), we thus constructed an 
equation of state of dense matter which describes a fluid 
of Skyrmions coupled to a dilaton field (associated with 
the glueball of QCD) and a vector meson field (coupled 
to the baryon number). Using our code, we constructed 
non-rotating, zero-temperature stable compact objects - 
we called Skyrmion Stars - and calculated their masses 
and radii to be 0.4 < M/Mq < 2.95 and 11.0 km < 
R < 15.3 km, respectively. In §5.1 in OB we compared 
Skyrmion stars - as defined in our model - to other 
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"Exotic" stars which also follow from solutions to an ef- 
fective non-linear field theory of strong forces (see also 
Heusler, Droz, & Straumann 1992 and references therein). 
Skyrmion stars in our picture are not boson/soliton stars 
where the soliton is a global structure over the scale of the 
star but rather form their constituent baryons as topolog- 
ical solitons using pions fields. 

In this paper, we compute models of rotating Skyrmion 
stars using the RNS code written and made publicly avail- 
able by N. Stergioulas (Stergioulas & Friedman 1995). 
RNS constructs models of rapidly rotating, relativistic, 
compact stars and assumes uniform rotation. The compu- 
tation solves for the hydrostatic and Einstein field equa- 
tions for uniformly rotating mass distributions, under the 
assumptions of stationarity, axial symmetry about the ro- 
tation axis, and reflection symmetry about the equatorial 
plane. The paper is presented as follows: In §0, we de- 
scribe the equations solved by the RNS code and then 
construct the corresponding rotating Skyrmion stars. The 
results are analyzed and discussed in §y. A comparative 
study of Skyrmion stars and neutron stars is done in ^ 
before concluding in §g|. 

2. Rapidly and rigidly rotating relativistic 
Skyrmion stars 

Here we present a brief outline of the equations solved for 
by the RNS code. We start by describing the space-time 
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around a rotating compact star in quasi-isotropic coordi- 
nates, as a generalization of Bardeen's metric (Bardeen 
1970): 



ds' 



^J+Pdt^ + e^'^ir^de^+dr^) 



(1) 



+ e 



7-p^2 



r''sm{ey{d(t)~ujdty 



(c = 1 = G). The metric potentials, 7, p, a, and the 
angular velocity of the stellar fluid to the relative local in- 
ertial frame (w) are all functions of the quasi-isotropic ra- 
dial coordinates (r) and the polar angle (9). The function 
(7+p)/2 is the relativistic generalization of the Newtonian 
gravitational potential; the time dilation factor between 
an observer moving with angular velocity co and an ob- 
server at infinity is cxp(l/2(7 + p)). We investigate uni- 
formly rotating perfect fluid stars with the energy momen- 
tum tensor given by: 



T'"' = (e + P)u''u'' + Pg"" 



(2) 



where e is the total energy density, P the pressure and u^ 
the unit time-like four velocity vector that satisfies 



-1 



(3) 



The proper velocity v of the matter, relative to the local 
Zero Angular Momentum Observer (ZAMO), is given in 
terms of the angular velocity fi = u^ /u^ of the fluid ele- 
ment (measured by a distant observer in an asymptotically 
flat space-time), by the following equation: 

v={VL-uj)rsui{e)e~P . (4) 

The four velocity (m'') of the matter can be written as 



g-(7+p)/2 
(1-1-2)1/2 



(l,0,0,f]) 



(5) 



Substitution of the above equation into Einstein field 
equations projected on to the frame of reference of a 
ZAMO yield three elliptic equations for the metric po- 
tentials p, 7 and Lo and two linear ordinary differen- 
tial equations for the metric potential a (Bardeen 1971; 
Butterworth, & Ipser 1976; Komatsu, Erigushi, & Hachisu 
1989). In general, the elliptic differential equations are 
converted to integral equations for the metric potentials 
using Green's function approach. 

From the relativistic equations of motion, the equa- 
tion of hydrostatic equilibrium for a barytropic fluid (the 
field equations are integrated by assuming various zero- 
temperature, barotropic EOS of the form e = e(p)) may 
be obtained as : 

h{P) -hp^-[jp + pp--f-p - ln(l - v^) 

+ A\cj~n,f], (6) 

where h{P) is termed as the specific enthalpy. Pp is the 
rescaled values for pressure, and hp is the specific enthalpy 
at the pole. 7p and pp are the values of the metric poten- 
tials at the pole, and fl = r^^. A is a rotation constant 
(Komatsu, Erigushi, & Hachisu 1989) and the subscripts. 



p, e and c stand for polar, equatorial and central, respec- 
tively. RNS solves the integral equation for p, 7 and lu, the 
ordinary differential equation (in 9) for the metric poten- 
tial a, together with Eq. (O), the hydrostatic equilibrium 
equations at the center [h(Pc) given] and at the equator 
[h{Pc) — 0], iteratively to obtain p, 7, a, to, the equato- 
rial coordinate radius (r^), angular velocity (fi), and the 
density (e) and pressure(P) profiles. 

The geodesic motion for circular orbits around rotat- 
ing relativistic stars is also presented in Bardeen (1970) 
and correspond to the motion of ZAMO yielding two pos- 
sible values for the velocity corresponding to corotating 
and counter- rotating orbits. The innermost stable circu- 
lar orbits are given by solving for V,rr — where V is 
the effective potential as written in Bardeen (1972); the 
comma followed by rr represent a second order partial 
derivative with respect to radius r. 

3. Numerical results 

3.1. Maximally rotating configurations 

The highest energy density used in all the computed ro- 
tating models is the value which gives the maximum mass 
non-rotating star as calculated in OB. We specify the 
equation of state, and the central energy density and the 
code computes models with increasing angular velocity 
until the star is spinning with the same angular velocity 
as a particle orbiting the star at its equatoin. 

Fig. nk and lb show the resulting stellar masses and 
radii as a function of the central density, respectively. The 
mass range is 0.4 < M/Mq < 3.45 while the correspond- 
ing radii (equatorial circumferential radius; [proper equa- 
torial circumference] /27r) are calculated to be 18.6 km < 
R < 23.0 km. In Fig. 1^, we show the resulting Mass- 
Radius plane. In general, for a given central density, rota- 
tion allows the masses and radii to increase by 30 % and 
40%, respectively, when compared with the non-rotating 
cases. The amount of baryonic mass in the outer region 
of the star (the crust region where p < pN is constructed 
using the EOS of Baym, Pethick, & Sutherland (1971); 
Pn is the nuclear saturation density) decreases drastically 
with rotation. Rotating Skyrmion stars crust constitute 
less than 5% of the total baryonic mass while it aver- 
ages 20% for the static configurations. The nuclear crust 
is further addressed in OB to which we refer the interested 
reader. 

The minimum spin period as a function of mass is 
shown in Fig. ||a and is calculated to be 0.8 ms < P < 
2.0 ms. The inertial frame dragging at the center of the 
star to the rotation rate is plotted in Fig. gb which shows a 
significant value throughout and reaches a maximum value 
of 0.8 for the most massive stars. It suggests that frame- 
dragging might have a significant effect on the structure 

We will sometimes refer to it as the Kepler frequency which 
is the frequency of a particle in a stable orbit at the circum- 
ference of a star. It is therefore also the frequency at which 
disruption of the star would occur. 
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Fig. 1. Maximally rotating configurations: a). 

Gravitational Mass versus central density (in units of the 
nuclear saturation density) for zero temperature Skyrmion 
stars in hydrostatic equilibrium. In this figure, the solid 
curves represent the rotating configurations while the 
dotted curves represent the non-rotating or static con- 
figurations, b). Radius versus central density, c). The 
Mass— Radius plane. 



Fig. 2. Maximally rotating configurations: a). 

Period versus gravitational mass for zero temperature 
Skyrmion stars in hydrostatic equilibrium, b). Ratio of 
the inertial frame dragging at the center of the star to the 
rotation rate versus gravitational mass. In this Figure and 
the subsequent ones, the vertical dotted line illustrates the 
mass of the maximum radius configuration. 



of rotating Skyrmion stars as P decreases since it is tak- 
ing place against the background of a radially dependent, 
frame-dragging frequency (see Weber & Glendenning 1992 
for example). 

It has been noted by Ha^nsel & Zdunik (1989) and 
Friedman & Ipser (1992) that the maximum spin rate for 
many neutron star EOS seems to be given by 



i £7] 



x( 



Afs., 



Mn 



.)l/2( 



R. 



10 km 



-3/2 



(7) 



where x has been quoted as either 7600 or 7700 s~^ . 
Ms, max. and i?s,max. are the total mass-energy and areal ra- 
dius of the maximum-mass static configuration for a given 
EOS. The newest fit for the proportionality constant can 
be found in Haensel, Salgado, & Bonazzola (1995) where 
X = 7600—7900 s^^ is given. If we assume that the mini- 
mum spin period for the Skyrmion star is 0.8 ms (that is 
the star with the maximum mass on the Kepler sequence 
has the maximum rotation rate) then the maximum an- 
gular velocity is 7854 rad s~'^ . Using the values for the 
maximum mass spherical star we get, 
^max = 7600-7900 X (2.95)1/2(1 44)-3/2 

= 7555-7852 s^^ . (8) 

This corresponds to less than few percent error which 
makes Eq. (0) an excellent tool for predicting the max- 
imum angular velocity even for our EOS. 



The centrifugal flattening of maximally rotating 
Skyrmion stars is obvious from Fig. pp. with eccentric- 
ities e = yj\ — [Rp/ReY {Rp and Re being the polar 
and equatorial radius of the configuration, respectively) 
ranging from 0.78 for the lightest stars up to 0.86 for the 
heaviest ones. The maximum eccentricity is reached at 
AI = 2.62Mq (the maximum radius configuration) after 
which e decreases, as expected. Fig. |p shows the Kinetic 
to gravitational energy ratio (T/|M^|) versus gravitational 
mass. Configurations with T/jT^/'l above the dashed line 
might be unstable to the bar mode. We come back to this 
in §3.2 where we discuss the conditions for the onset of 
non-axisymmetric modes and to what extent they could 
be used as a stronger constraint than the maximum rota- 
tional frequency. 

In Fig. H, the moment of inertia (Fig. 0a) and the angu- 
lar momentum (Fig. 0b) are plotted. The sudden decrease 
in the moment of inertia for M > S.SM© is the conse- 
quence of the mass reaching a "plateau" (at M = 3.3Mq) 
while the radius keeps decreasing (see Fig. |l|a and Fig. |l|b) . 
The combined effects of / decreasing with Q increasing 
(Fig. ya) explains then why the angular momentum stays 
constant for M > 3.3Mq - using the simple dimensional 
analysis J = I ft. Finally the height from surface of the 
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Fig. 3. Maximally rotating configurations: a). 

Eccentricity versus gravitational mass for zero temper- 
ature Skyrmion stars in hydrostatic equilibrium, b). 
Kinetic to gravitational energy ratio versus gravitational 
mass. The dashed line is based on the empirical formula 
for the onset of the bar mode instability as described by 
Eq.p. Configurations with T/|VF| above the dashed line 
might be unstable to the bar mode. 



last stable co-rotating orbit (h^ = Torb. — R)q is shown in 
Fig. 0c. It can be seen from this figure that in general the 
stable orbits exists up to the surface of the star. The large 
radii calculated makes it very difficult for the stable orbit 
to be outside the star (when compared to 6GM/c^; the 
inner most stable Kepler orbit around the non-rotating 
case). For the most massive configurations (M > 3.2Mq), 
the boundary layer (the separation between the surface of 
the neutron star and its innermost stable orbit) can be as 
as high as 1.5 km for the maximum value. Lowering the 
spin frequency will eventually relax the above conclusions 
(see §3.3). 

3.2. Bar instability 

Skyrnrion stars are subject to rotational constraints linked 
to gravitational- wave instabilities which make all rotating, 
perfect fluid equilibria unstable to modes with angular 
dependence expiimcj)) for sufficiently large m (Friedman 
& Schutz 1978; Lindblom 1984). These instabihties allow 
the star to convert its rotational energy to gravitational 
energy waves (Glendenning 1997). 



Fig. 4. Maximally rotating configurations: a). 

Moment of inertia versus gravitational mass for zero tem- 
perature Skyrmion stars in hydrostatic equilibrium, b). 
Angular momentum (Dimensionless ratio of angular mo- 
mentum cJ/GMq) versus gravitational mass. c). Height 
from surface of last stable co-rotating orbit in equatorial 
plane versus gravitational mass. 



Fairly recent general relativistic calculations have 
shown that the gravitational-radiation-driven bar-mode 
(/-mode) instability sets in at values of T/|M^| much lower 
than in the Newtonian limit. For the case of realistic EOS, 
Morsink, Stergioulas, & Blattnig (1999) found the follow- 
ing empirical formula for the onset of the bar mode insta- 
bihty (the case of polytropes can be found in Stergioulas 
& Friedman 1998): 



T/\W\ = 0.115 - 0.048M/MX'L 



(9) 



where M^^^ is the maximum mass star for the non- 
rotating case. Assuming that this formula holds for our 
EOS, we find that all maximally rotating stars with RI > 
Mq might be prone to such an instability. 

The real question, however, is whether the /-mode 
instability is actually going to limit the spin period of 
Skyrmion stars. With the onset of the m — 2 (/-mode), 
the r-modeR will soon become unstable (Friedman & 
Morsink 1998) and spin down the star to even slower rate. 
As a result, it is not expected that the /-mode will limit 
the spin-rate of any stars (Andcrsson, Kokkotas, & Schutz 



For the cases where rorb. is non-existent or inside the star, 
Torb. is taken to be the Keplerian orbit radius at the surface of 
the star. 



■^ The r-mode instability is a member of the class of 
gravitational-radiation-driven instabilities and are generic in 
rotating neutron stars (Friedman & Morsink 1998). That is, 
every r-mode is in principle unstable in every rotating star, in 
the absence of viscosity. 
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1999), although it may become unstable and emit gravi- 
tational radiation in some cases. 



3.2.1. Skyrmion fluid viscosity 

Realistic compact stars are viscous, and the presence 
of viscosity will shift the onset of instability and if 
large enough will damp out the instability (Andersson, 
Kokkotas, & Schutz 1999). If the instability is driven by 
viscous dissipation, it was shown in Bonazzola, Frieben, & 
Gourgoulhon (1995) that for cold uniformly rotating NS 
only very stiff equations of state arc likely to allow for 
spontaneous symmetry breaking. If we consider the stiff- 
ness of our EOS (when compared to standard EOS; see 
Figure 2 in OB) then instability to a bar mode is a very 
likely plausibility in Skyrmion stars. So how viscous is a 
Skyrmion fluid and what is its temperature dependency? 

The Skyrmion fluid can be looked at as made of 
fermionic soliton objects and in principle one should be 
able to calculate its physical parameters, among others its 
viscosity. The complication is linked to the fact that un- 
like a pure fermionic fluid, the structure of the Skyrmion 
fluid is highly non-linear (a consequence of the Skyrmion 
having structure). This makes the calculations not trivial 
and beyond the scope of this paper. There might still be 
the possibility that such a fluid is viscous enough that it 
will of course damp these instabilities altogether so that 
limits on rotation is set by the Kepler frequency - a notion 
to be confirmed. 



0.9 
0.8 
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Fig. 5. Different spin frequencies: a). The ISCO fre- 
quency versus mass for Skyrmion stars rotating at 290, 
360 and 580 Hz respectively. The dotted line represents 
the static models, visco = 2198 Hz(M0/M). b). Height 
from the star surface of the innermost stable circular or- 
bit (ISCO) versus mass for Skyrmion stars rotating at 290, 
360 and 580 Hz respectively (see Table 1). 



3.3. Different spin frequencies: Skyrmion stars and 
QPOs 

The discovery of kilohertz Quasi-Periodic Oscillations 
(QPOs) in low mass X-ray binaries (LMXBs) with the 
Rossi X-Ray Timing Explorer has stimulated extensive 
studies of these sources (van der Klis 1998). One thing that 
has been suggested in the literature is that the highest fre- 
quency QPO observed could correspond to the orbital fre- 
quency at the inner stable circular orbit (ISCO) (Miller, 
Lamb, & Psaltis 1998). The signature of the marginally 
stable orbit is a saturation in QPO frequency (assumed to 
track inner disk radius) versus mass accretion. Such a sat- 
uration where the frequency becomes independent of the 
chosen mass accretion rate indicator has been reported in 
4U 1820-30 (Zhang et al. 1998; Kaaret et al. 1999). 

In Fig. 1^ we plotted ISCO frequency (for the prograde 
case only) versus Mass for Skyrmion star models rotating 
at frequencies 290, 360 and 580 Hz respectively. These val- 
ues allow us to cover the range of frequencies observed in 
LMXBs. Fig. ^3 shows the corresponding gaps; that is the 
height between the star surface and the ISCO (prograde 
case). In general, the Skyrmion star mass range with an 
existing gap is calculated to be. 



while the corresponding radii and visco are 



14.0 km < i? < 20.0 km 



(11) 



0.8 kHz < j/isco < 1-3 kHz 



(12) 



respectively. Our numbers can thus account for the high 
masses (up to 1.78 ± O.23M0; Orosz, & Kuulkers 1999) 
determined from LMXBs (see also Miller, Lamb, & Cook 
1998). General features of the models discussed in this 
section are shown in Table 1. A key to all the the tables 
in this paper is as follows: 



1.8 Mq < M < 3.0 M, 







(10) 
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Table 1. Skyrmion stars: Different spin frequencies (vs) 
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0.00 


0.31 


2.36 
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0.016 


0.31 


5.0x10" 


2.28 


2.59 


16.47 


1.93 


6.26 


0.28 


5.09 


1.06 


0.011 


0.49 


6.9x10" 


2.60 


3.02 


16.13 


4.89 


9.31 


0.24 


5.83 


1.21 


0.009 


0.59 
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: 360.0 H 
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S.OxlO"-* 


1.38 


1.48 


15.93 


0.00 


0.00 


0.39 


2.48 


0.64 


0.025 


0.31 


5.0x10" 


2.30 


2.62 


16.62 


1.69 


7.08 


0.34 


5.24 


1.35 


0.017 


0.50 


6.9x10" 


2.62 


3.04 


16.24 


4.53 


10.04 


0.28 


5.96 


1.53 


0.014 


0.59 


1.0x10^^ 


2.80 


3.29 


15.26 


7.06 


12.18 


0.24 


5.81 


1.50 


0.010 


0.69 
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3.0x10"'' 


1.53 


1.66 


17.43 


0.00 


3.33 


0.64 


3.25 


1.35 


0.069 


0.34 


5.0x10" 


2.43 


2.76 


17.42 


1.31 


10.17 


0.53 


6.05 


2.51 


0.048 


0.52 


6.9x10" 


2.72 


3.16 


16.78 


3.58 


12.72 


0.47 


6.58 


2.73 


0.038 


0.61 


1.0x10^^ 


2.86 


3.36 


15.57 


5.85 


14.29 


0.39 


6.18 


2.56 


0.028 


0.70 



Cc 

Mg 
Mo 
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'^orb 

h+ 

h- 

e 
I 

J 

T/\W\ 

Zp 
Zf 



Star's angular velocity (10* Hz) 

Star's spin frequency, v = f7/27r (Hz) 

Central density (g cm^'^) 

Gravitational mass (Mg/Mq) 

Baryonic mass (Mq/Mq) 

Radius of star, measured at equator (km) 



bers remain to be confirmed from observations and furthet 
modelling of the system. 

3.3.2. The case of 4U 0614+09 

The highest known QPO has a frequency of 1329 Hz (van 



Radms of the mnermost stable orbit (km) „, , , , „„„„n mi r n.i t 

-TT • , i r ii i r r • / _L ibtraaten et al. 2000). Ihe frequency of the correspondmg 

Height from the star surface of innermost statne , ,-,^ „^, , , „„^ . , ,. 7 , , ^^r, -r-, r^ i 

prograde circular orbit (km) 

Height from the star surface of innermost stal; 

retrograde circular orbit (km) 

Eccentricity 



Moment of Inertia (10''^ g cm^) 
Angular momentum (cJ/GMq) 
Rotational energy to the gravitational energy 
Ratio of the inertial frame dragging at the 
center of the star to the rotation rate 
Polar redshift 
Forward redshift 
Backward redshift 



To first order in the rotation rate of the star {i/g < 400 
Hz), the orbital frequency in the marginally stable orbit 
is Pisco ^ 2198 Hz(M0/M)(l + 0.748j(Mq/M)2) (see 
Miller, Lamb, & Cook 1998, for example). Consider, as 
an example, a moderately slowly rotating Skyrmion star, 
1^ = 290 Hz, with mass M — 2.6Mq. Using the value of 
j given in Table 1 we find visco ~ 9589 Hz while the 
value in the non-rotating case is 8454 Hz (depicted by the 
dotted line in Fig. ||a). 

3.3.1. The case of 4U 1820-30 

We first apply our model to 4U 1820-30 which seems to 
present the strongest experimental evidence for the exis- 
tence of the marginally stable orbit (Zhang et al. 1998; 
Kaaret et al. 1999). For t'isco — 1060 Hz as the estimated 
frequency at the ISCO, and the spin frequency to near 
300 Hz, Fig. ||a suggests a mass 2.35M0. The radius is 
then calculated to be of the order of 16 km. These num- 



jle 

sfar (4U 0614-^09) is beheved to be near 300 Hz. Simple 

itimates (van Straaten et al. 2000, §4.4; van der Klis 2000, 
k6) show the difficulty for extremely stiff EOS such as 
ours to account for the derived range in mass and radius 
(M < 1.9M0 and R < 15.2 km). We have noted that for 
our EOS visco < 1-3 kHz and i? > 14 km. Recent model- 
ing put even more constraints by suggesting a more com- 
pact 4U 0614+09 with a radius less than 10 km (Titarchuk 
& Oshcrovich 2000). 

It is not trivial to account for such numbers in our 
model unless the 1329 Hz frequency does not actually 
correspond to the ISCO. Indeed, unlike in the 4U 1820- 
30 source, there is no obvious saturation of the kilo-hertz 
QPO with respect to mass accretion indicator in the 4U 
0614-^09 case (see Fig. 4 in van Straaten et al. 2000). 



4. A comparative study of Skyrmion stars and 
Neutron stars 

In this chapter we compare Skyrmion stars to the more 
traditional neutron stars. We chose three recent/modern 
EOSs: the first is described in Baldo, Bombaci, & Burgio 
(1997) which uses the Argonne vl4 (Avl4)) (Wiringa, 
Smith, & Ainsworth 1984) two-body nuclear force; the sec- 
ond uses the Paris two-body nuclear force (Lacombe et al. 
1980) , implemented in both cases by the Urbana three- 
body force (TBF) (Carlson, Pandharipande, & Wiringa 
1983; Schiavilla, Pandharipande, & Wiringa 1986) while 
the third one is that developed in Prakash et al. (1997) 
using a generalized Skyrme-like EOS. As in Datta, 
Thampan, & Bombaci (1998), we refer to these EOSs as 
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Fig. 6. The M — R relation for non-rotating Skyrmion 
stars (OB) as compared to theoretical models of non- 
rotating neutron stars (UU, BBBl, BBB2, BPAL12, Hyp, 
and K~^) and strange stars (SSI and SS2). The data for 
the neutron stars and strange stars was kindly provided to 
us by Li et al. (1999). The Schwarzschild radius {2GM/c^) 
is shown as a dotted line. Inside the triangle is the allowed 
range of M and R for 4U 1636-53 as modeled in Nath, 
Strohmayer, & Swank (2001) using fits to X-ray bursts. 

BBBl (Avl4+TBF), BBB2 (Paris+TBF) and BPAL32, 
respectively. 

Table 2. Maximum-mass non-rotating models 



EOS 


Ec 


Mg 


Mo 


R 


rorb 


BBBl 


3.09x10"'' 


1.788 


2.082 


9.646 


15.845 


BBB2 


3.12x10^^ 


1.917 


2.261 


9.519 


16.984 


BPAL32 


2.67x10^^ 


1.947 


2.263 


10.509 


17.254 


OB 


1.22x10"^ 


2.797 


3.299 


14.565 


24.743 



Table 2 summarizes the non-rotating compact star 
structure parameters for the EOS models BBBl, BBB2, 
BPAL32 and ours (OB). The values listed correspond to 
the maximum stable mass configuration. Note that for this 
table and the rest of the tables, we consider only the co- 
rotating case (the cases where rorb. is non-existent or in- 
side the star, rorb. is taken to be the Keplerian orbit radius 
at the surface of the star). The maximum mass is usually 
an indicator of the softness/stiffness of the EOS and its 
values as listed in Table 2 confirms our previous state- 
ments that our EOS is stiffer than the standard ones. The 
consequences of such a stiffness is also illustrated in Table 
3 where we list the quantities corresponding to the max- 
imum gravitational mass configurations. It can be seen 



that the gravitational mass of the maximum stable ro- 
tating Skyrmion star has a value of ~ 3.5Mq, while it 
is ~ 2AMq for the configurations constructed with the 
other EOS. That is, with our EOS we do not need to go 
to break up speeds to account for the ~ 2.0Mo mass pre- 
dicted from analysis of LMXB observational data (Zhang 
et al. 1996). The values of these same quantities are listed 
in Table 4 for the maximum angular momentum models 
leading to similar conclusions. 

4.1. 4U 1636-53 as a Skyrmion star candidate? 

In Fig. ra wc compare the M — R relation for Skyrmion 
stars (OB) to the theoretical M — R curve obtained us- 
ing six recent realistic models for the EOS (UU, BBBl, 
BBB2, BPAL12, Hyp, and K"!). The sohd curves labeled 
SSI and SS2 are for strange stars (the data was kindly 
provided to us by Li et al. 1999). The triangle depicts 
the mass-radius constraint from fits to X-ray bursts in 
4U 1636-53. Inside the triangle is the allowed range of 
M and R which satisfies the compactness constraints as 
modeled in Nath, Strohmayer, & Swank (2001; see their 
Figure 4), and clearly favoring stiffer EOSs. Our modeled 
stars (OB) cross the triangle suggestive of 4U 1636-53 as 
a plausible Skyrmion star candidate. However, one should 
keep in mind the fact that modern EOSs can be modified 
as to also cross the triangle (Heiselberg & Hjorth-Jensen 
1999). 

5. Conclusion 

With the RNS code, we constructed numerical models 
of rotating Skyrmion stars for a newly derived EOS of 
dense matter based on a fluid of Skyrmions. We calcu- 
lated their masses and radii to be 0.4 < M/Mq < 3.45 
and 13.0 km < R < 23.0 km, respectively. The spin period 
of the maximally rotating configurations are calculated to 
be 0.8 ms < P < 2.0 ms. Owing to the stiffness of the 
Skyrmion fluid, Skyrmions stars are found to be on aver- 
age heavier, to show higher equatorial radii, and to rotate 
slower than theoretical models of neutron stars based on 
modern realistic EOS. 

Skyrmion stars might have little bearing on reality 
since there is still no guarantee that our EOS (OB) is a 
plausible representation of matter above nuclear densities. 
However, our model so far has been successful in reproduc- 
ing the basic features of compact objects. An interesting 
consequence of our model, with its plausible applications 
to QPO systems, is the fact that massive Skyrmion stars 
(1.8 < M/Mq < 3.0) can possess gaps with orbital fre- 
quencies in the kHz range (0.8 kHz < t'isco < 1-3 kHz). 
These points suggest that the model warrants further 
study. 
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Table 3. Maximum-mass rotating models 



EOS 



Mg Mo 



R 



'orb 



BBBl 2.56x10''' 

BBB2 2.82x10^'^ 

BPAL32 2.27x10^'^ 

OB 1.03x10^" 



2.135 
2.272 
2.300 
3.445 



2.471 
2.653 
2.657 
4.034 



13.129 13.490 

12.519 13.550 

14.276 14.611 

19.890 21.089 



0.703 
0.687 
0.699 
0.837 



Table 4. Maximum angular momentum models 



EOS 



Ma Mo 



R 



'orb 



BBBl 

BBB2 

BPAL32 

OB 



2.44x10'" 
2.82x10^^^ 
2.14x10^^^ 
9.00x10" 



2.133 
2.272 
2.299 
3.435 



2.468 
2.653 
2.655 
4.017 



13.264 13.558 

12.519 13.550 

14.481 14.713 

20.549 21.498 



0.706 
0.687 
0.702 
0.840 



n. 



1.095 
1.203 
1.001 
0.749 



I 


J 


T/W 


UJc/i^c 


z. 


Zf 


Zb 


2.428 


3.019 


0.120 


0.764 


0.690 


-0.330 


1.975 


2.539 


3.469 


0.123 


0.825 


0.849 


-0.349 


2.483 


3.005 


3.416 


0.113 


0.771 


0.679 


-0.328 


1.933 


9.982 


8.507 


0.144 


0.777 


0.777 


-0.342 


2.287 



n. 



1.079 
1.203 
0.981 
0.715 



/ 


J 


T/W 


UJc/f^c 


z. 


Zf 


Zb 


2.465 


3.021 


0.120 


0.756 


0.677 


-0.328 


1.935 


2.539 


3.469 


0.123 


0.825 


0.849 


-0.349 


2.483 


3.070 


3.419 


0.114 


0.760 


0.661 


-0.326 


1.881 


10.52 


8.560 


0.146 


0.751 


0.728 


-0.355 


2.133 
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